Relations between three formalisms 
for irrotational binary neutron stars in general relativity 
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Various formalisms proposed recently for irrotational binary systems in general relativity are 
compared and explicit relations between them are exhibited. It is notably shown that the formalisms 
of (i) Teukolsky, (ii) Shibata and (iii) Bonazzola et al. (as corrected by Asada) are equivalent, i.e. 
yield exactly the same solution, although the former two are simpler than the latter one. 
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I. INTRODUCTION 



Inspiraling neutron star binaries are expected to be one of the strongest source of gravitational radiation for the 
interferometric detectors GEO600, LIGO, TAMA and VIRGO currently under construction . These systems are 
therefore subject to numerous theoretical studies. Among them are fully relativistic hydrodynamical treatments, 
pioneered by the works of Wilson et al. jlJJ^I and Oohara and Nakamura j|. The most recent numerical calculations, 
| those of Baumgarte et al. and Marronetti et al. ||, rely on the approximations of (i) quasiequilibrium state 
■ and (ii) synchronized binaries. Whereas the first approximation is well justified up to the innermost stable orbit, the 
second one does not correspond to physical situations, since the gravitational-radiation driven evolution is too rapid 
for the viscous forces to synchronize the spin of each neutron star with the orbital motion |7j|| — as they do for 
ordinary stellar binaries. Rather, the fluid velocity circulation (with respect to some inertial frame) is conserved in 
these systems. Provided that the initial spins are not in the millisecond regime, this means that close configurations 
are well approximated by irrotational (i.e. with zero vorticity) states. 

The first relativistic formulation for quasiequilibrium irrotational binaries has been given by Bonazzola, Gourgoulhon 
and Marck H (hereafter BGM). Their method is based on one aspect of irrotational motion, namely the counter- 
rotation (as measured in the co-orbiting frame) of the fluid with respect to the orbital motion. Then Teukolsky |Tc| ] 
and Shibata jllj gave two other formulations based on the very definition of irrotationality, which implies that the 
specific enthalpy times the fluid 4- velocity is the gradient of some scalar field (potential flow) . The formulation 
... of Teukolsky is fully four-dimensional, whereas that of Shibata is three-dimensional. 

The aim of the present note is to clarify the relations between these three formalisms and in particular to answer 
to the question: do these rather different formalisms lead to the same solution ? Of course, at the Newtonian limit, 
they do. However it is not obvious that this still holds in the relativistic regime. For instance, the first integral of 
motion — the fundamental equation which is used to compute the matter distribution in numerical procedures — is 
written by BGM as (Eq. (66) of Ref. @) 

In h + - ln(N 2 - BiB 1 ) + In T = const , (1) 

whereas for Teukolsky it reads (Eq. (54) of Ref. [jToj] ) 

N[h 2 + Dill) L>V] 1/2 + B i D i ip = const =: C (2) 
and for Shibata (Eq. (2.18) of Ref. 0) 

^ + S l D^ = const . (3) 
A 

In these equations, h is the fluid specific enthalpy, N is the lapse function corresponding to the 3+1 slicing of spacetime 
by spacelike hypersurfaces T, t , B l is the shift vector of rotating coordinates^], T is the Lorentz factor between the 



1 Greek (resp. Latin) indices run from to 3 (resp. 1 to 3). Unless explicitly mentioned, the notations of BGM are 
systematically used. In particular, the shift vector B l is the negative of B % defined by Teukolsky JIG]. 
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fluid observer and the co-orbiting observer, ip is the scalar potential from which the fluid 4-velocity can be derived 
(irrotational condition) , D t is the covariant derivative on the hypersurface (related to the full spacetime covariant 
derivative V Q by Eq. (33) of Ref. pOfl ), and A and S l define the decomposition of the fluid velocity between a part 
along the helicoidal Killing vector of the quasiequilibrium assumption || and a part parallel to E t 0. Needless to say, 
Eqs. (0)-(||) differs from each other quite substantially, at least at first glance. However, we shall see that they do 
represent the same equation. More generally, the answer to the question raised above is t/es, i.e. BGM, Teukolsky 
and Shibata formalisms are equivalent. To establish this, we will first consider Teukolsky formulation (Sect. |l|) and 
make the connection of the two other formulations with it (Sect. |l| and IV). 



II. TEUKOLSKY FORMALISM 

Let us briefly recall and discuss Teukolsky formalism |Io| . The neutron star matter is very well represented by a 
perfect fluid, whose stress-energy tensor writes T a p = (e + p)u a u/3 + pg a p, e being the fluid proper energy density, p 
the fluid pressure, u a the fluid 4-velocity and g a p the spacetime metric. The zero temperature approximation is fully 
justified for neutron star matter. In this case the fundamental energy-momentum conservation relation V M T MQ = 
can be shown to be equivalent to the two equations 

u t "V fl (hu a ) + W a h = , (4) 

V M (nu") = 0, (5) 

where n is the baryon numerical density and h is the specific enthalpy: h = (e+p)/ (mbfi), being the baryon mass. 
Note that thanks to the First Law of Thermodynamics at zero temperature, h is also the baryon chemical potential 
(divided by m^), so that the Gibbs-Duhem relation writes S/ a p = mbuV a ft, relation which has been used to derive 
Eq. i). 

The vorticity, or rotation 2-form, of the fluid is defined as |13| | 

^=P^P/V [ ^ ] , (6) 

where P a p := g a p + u a up is the projection tensor on the fluid rest-frame hyperplanes. Note that the definition (||) 
differs from that of Teukolsky (Eq. (19) in Ref. |l(|), which includes a factor h in front of u a . However, it is easy to 
see that VL^^ ^ 7 — 2htt a p, so that an irrotational flow can be defined either by Q a p — or by r2^ ukolsky = 0. By 
means of Eq. ([|) , the vorticity can be re- written as 

tt a/ 3 = ~V[p(hu a ]) , (7) 

so that the irrotational condition simply means that the 1-form h u a is closed. Hence there exists (at least locally) a 
scalar field tp such that 

u a = , (8) 

i.e. the flow is potential The normalization relation u^u^ = — 1 implies then 

h 2 = -V^V^. (9) 

Thanks to this latter equation, the momentum conservation equation (^|) is automatically satisfied, so that the problem 
reduces to finding a solution of Eq. (|^), which becomes 

V^Vip + Vif) V M \n{n/h) = . (10) 

The quasiequilibrium assumption is implemented by requiring the existence of a helicoidal Killing vector l a ||. 
Teukolsky has then shown that Eq. (g) is equivalent to Eq. (||) and that Eq. ( |io| ) becomes (Eq. (50) of Ref. []l0|| ) 



A and 5* 1 are respectively denoted u° and V % by Shibata [ jnj . 
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A-DV + £>ty A In ^ j + ^2 (C - A* A^) A In (^) - JVif) - Di{B* Drf) = , (11) 

where K is the trace of the extrinsic curvature tensor of the S t hypersurfaces and C is the same constant as in 
Eq. (H). C is in fact defined by Eq. (39) of Ref. [jTo| : C = —l^V^ip. The three-dimensional equation ( |Tl|) has to be 
solved in -0 to get a solution to the problem. 



III. LINKS BETWEEN SHIBATA AND TEUKOLSKY FORMALISMS 

The starting point of Shibata formulation |ll]] is the following decomposition of the fluid 4- velocity (see footnote^): 

u a = X(l a + S a ) with n^S" = , (12) 

where n a is the unit future-directed normal vector to S t . This relation is inserted in the momentum conservation 
equation and the result projected onto £ t to get^| (Eq. (2.13) of Ref. |ll]] with the helicoidal symmetry taken into 
account) 

A (~ + Pl S l ) + Si(D jPi ~ Dfpj) = , (13) 



where 



P a := h a (hu») , (14) 



h a /3 '■= 9ap + n a np being the projection tensor onto St (p l is denoted by Shibata). Shibata defines then the 
irrotational condition by requiring 

Pi = Di(j> , (15) 



where is a scalar field. Equation (13) reduces then to the first integral (J3j). By a somewhat lengthy calculation, 
Shibata shows that ( |l5| ) is equivalent to f2 a/ g = 0. As discussed in Sect, [fi|, this means that the definitions of 
irrotationality by Shibata and Teukolsky coincide. It must then be possible to exhibit a correspondence between the 
two formulations. 

First, inserting Eq. (||) into Eq. ( |l4]) and using Eq. Jl|) leads immediately to <fi = ijj. Then, the expression 
l a = Nn a - B a (Eq. (7) of Ref. f|) combined with Eq. (|lj) leads to 

S a = \h a u" + B a = 1 - T -D a ^ + B a . (16) 
Besides, the normalization relation u^u fL = — 1 gives 

1/1 



X= N { 1+ V D ^ D ^) ■ 

Inserting Eqs. ([lii]) and ([l7|) in Shibata's form of the first integral of motion [Eq. (^) above] gives Teukolsky 's version 
of it [Eq. (§) above]. 

The equation for determining tp given by Shibata (Eq. (2.22) of Ref. Ell) is also at first glance quite different from 
the equation given by Teukolsky (Eq. ([ll]) above), since it reads 

Di (jV^Dty) + D i (NnXB i ) = . (18) 

However, Eqs. ( p~7| ) and (|^) yield the following expression of A: 

X = j^ K (C-B i D i ^) , (19) 

which, once reported in Eq. (|l8|), gives exactly Eq. (fir]). 



3 Note that Eqs. (2.10) and (2.11) of Ref. |n) are not correct: their right-hand-side must be supplemented by respectively 
—hn^u^KijV 3 and +hn l± u' J ' KijV 3 (using Shibata's notation for V 3 ). Fortunately, these corrections cancel each other when 
summing Eqs. (2.10) and (2.11), so that Shibata's final result (Eq. (2.13) of Ref. Jllf ) is correct. 
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IV. LINKS BETWEEN BGM AND TEUKOLSKY FORMALISMS 



The BGM formulation H is based on the fluid 3-velocity with respect to the co-orbiting frame. This later frame 
is that of the observer whose worldlines are the trajectories of the Killing vector l a , i.e. whose 4- velocity is v a = 
e~*/ Q , with $ = l/21n(iV 2 — BiB 1 ). The 3-velocity V a with respect to that observer is defined by the orthogonal 
decomposition 

u a = T(V a + v a ) with v^V 1 = . (20) 

The momentum conservation equation (||) is then re-expressed as an equation for V a , which is a relativistic gener- 
alization of Euler equation in a rotating frame. The basic idea of BGM is to define a counter-rotating motion as a 
motion for which the Coriolis term is canceled by the (V A V) /\ V term coming from the advection term V ■ W, 
where the wedge product is defined in the co-orbiting observer rest frame: (X A Y) a = v a ( aa llv X ll Y 1 ' . It has been 
however noted by Asada |l4j that the condition given by BGM is not enough to fully determine the velocity field. In 
fact the correct condition is obtained by setting F — in Eq. (50) of BGM ||, so that this latter becomes 

(V A e~^V) a — —2e~^u a , (21) 

which is identical to Eq. (4.12) of Asada |Q. In this equation u> a is the rotation vector of the co-orbiting observer: 
2uj a := v a t aa ^ v V ' uVw If the counter-rotation condition ( pl|) is satisfied, the Euler equation reduces then simply to 
the first integral ([l]). The baryon number conservation equation (g) gives an additional equation, for the divergence 
oiV a : 



+ ln(nr) = . (22) 

The BGM procedure is to perform a 3+1 decomposition of Eq s. (pl| )-(p2[) and to introduce scalar and vector potentials 
so that Eq. ( |22| ) gives a scalar Poisson-type equation and Eq. ([2l|) gives a vector Poisson equation. This contrasts with 
Teukolsky or Shibata formalisms, which require to solve only one scalar Poisson-type equation (Eq. ( pi] ) or Eq. (p"8|)). 
It has been however noted by Asada |L4j that the counter-rotation condition ( pl| ) is equivalent to the irrotational 
condition f2 Q ^ = 0. Hence there must exist a re-interpretation of the BGM formulation in terms of Teukolsky's one. 
In fact, if the irrotational condition (||) holds, one can see easily, by means of the relation T — —v^u 11 , that 

r = ;&- < 23 > 

where C is the same constant as in Eq. (|J). Taking the logarithm of expression (^3|) gives directly the first integral 

©■ ^ rn r-v 

Combining Eqs. @, (§) and (||) yields 

<E> 

V a = e —^ a i) - v a , (24) 
from which one can derive the following expression for the part of V a parallel to Et (cf. Sect. V.B of BGM): 

W i = — Dty + e-*B* . (25) 
o 

It is immediate to verify that the 3-velocity given by Eq. (|24|) obeys to the counter-rotation conditio n (|2l|) . Note 
that this way of proceeding is more straightforward than the converse one: it took a full appendix in [|14[to verify 
that counter-rotation implies irrotational flow. 



V. CONCLUSION 



We have shown explicit relations between the formalisms introduced by BGM, Teukolsky and Shibata to treat 
irrotational neutron star binaries. All these formalisms are equivalent in the sense that they lead to exactly the 
same solution. Teukolsky's formulation jTo| is the simplest one because it is directly based on the four-dimensional 
definition of irrotational motion. Shibata's one [Tl| ] contains some complications induced by its three-dimensional 
character, but both formulations result in the same type of scalar Poisson-like equation governing the fluid motion 
(Eqs. (11) and ( |l8| ) above). On the contrary, BGM formulation ||] (corrected by Asada |Q) is more complicated 
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since it requires the resolution of an additional vector Poisson equation. Therefore, for numerical studies, Teukolsky 
or Shibata procedure should be preferred. 
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